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Abstract 

We investigate the effects of the nonminimal coupling between the 
scalar field dark energy (quintessence) and the dark matter on the two- 
point correlation function. It is well known that this coupling shifts 
the turnover scale as well as suppresses the amplitude of the matter 
power spectrum. However, these effects are too small to be observed 
when we limit the coupling strength to be consistent with observations. 
Since the coupling of quintessence to baryons is strongly constrained, 
species dependent coupling may arise. This results in a baryon bias 
that is different from unity. Thus, we look over the correlation function 
in this coupled model. We are able to observe the enhancement of 
the baryon acoustic oscillation (BAO) peak due to the increasing bias 
factor of baryon from this species dependent coupling. In order to 
avoid the damping effect of the BAO signature in the matter power 
spectrum due to nonlinear clustering, we consider the coupling effect 
on the BAO bump in the linear regime. This provides an alternative 
method to constrain the coupling of dark energy to dark matter. 

Due to the strong constraint on the coupling of the scalar field dark 
energy (quintessence) to baryons from the local gravity, we investigate the 
efi^ect of the species dependent coupling [1, 2] by considering a model in 



which the quintessence Q is only coupled to the cold dark matter (CDM) . We 
assume a Yukawa-type coupling, rric = e^'^'^rn*, where m* is the bare mass 
of the CDM [3, 4]. This specific choice of coupling requires that the present 
value of the scalar field vanishes in order to satisfy rric = m* at present. 
Then we are able to write the general action including this interaction as 



S = -J d^xy^}^ — [-R + di'Qdi.Q] + V{Q) -Cc-Cr-Cb^, (1) 

where M = I/VSttG is the reduced Planck mass, V{Q) is the potential 
of Q, and C,- and Ci = —miS{x — Xi{t))^^ g^^.x'^x'^ / g {i = c,b) denote 
the Lagrangian of radiation, CDM, and baryons respectively. We adopt 
V{Q) = Vq exp(A(5^/2) with A = 5 in the following [4, 5]. However, the main 
conclusions are independent of the form of the scalar field potential (see be- 
low). Due to the coupling, the scalings of the CDM and the quintessence 
energy densities are changed respectively to [5] 

Pc{a) = pla-^+' where eln(a) = nc[Q(a)-Q(l)] , (2) 

p'q = -3Hil+u'^)pQ where a;gf = a;Q + ^-^Q', (3) 

^ H PQ 

where Ti = {da/dri)/a, loq is the equation of state (eos) of the quintessential 
dark energy (DE) , p° denotes the present value of the CDM energy density, 
the present value of scale factor ao = 1, and primes mean the differentiation 
with respect to the conformal time r]. Generally, the sign of e depends 
on both the model and the form of the coupling. The linear perturbation 
equations for the CDM and the scalar field Q in the synchronous gauge are 
[6] 

5c = -ec-^h' + ncSQ' , (4) 

e'^ = -Hec + nc(k^SQ-Q'ec) , (5) 

SQ" + 2HSQ' + eSQ + ^^SQ = -^h'Q'-^p,nJ,, (6) 

where k is the wave-number, h is the metric perturbation, Sc = Spc/pc, 
and 9c is the gradient of the CDM velocity fiow. Also from the perturbed 
Einstein equations, we obtain 

^{h"+Hh') = -^[2Spr+Spc+{M4)5p,+^Q'6Q'-2^dQ] , (7) 
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where C(, is the sound speed of baryons. Note that we will adopt the adiabatic 
initial conditions and thus k'^5c term is absent in Eq. (5). 

The coupling strength ric is commonly constrained through the compari- 
son with the observed cosmic microwave background anisotropy and matter 
power spectrum. In Ref. [5], we found ric < 0.01. Even though the actual 
value of the upper limit depends on the form of the quintessence potential 
and that of the coupling, the obtained limits for other potentials and cou- 
plings are of the same order as shown in Ref. [7]. Furthermore, Eqs. (4) 
and (5) , which describe the evolutions of the CDM density and velocity field 
respectively, imply that the influence of the quintessence field are dwarfed 
by the background evolution, %. Thus, as long as we have the latc-time 
dominated quintessence model, the evolution behaviors of Eqs. (4) and (5) 
are quite similar and almost independent of the form of potential. However, 
we have also checked the early dark energy model, for example, with the po- 
tential given by V{Q) = Vbcosh(A(5), in which the dark energy component 
is not negligible at early times [8]. In this model, % is quite different from 
the quintessence model with the exponential potential and produce a quite 
different behavior of 5c- Here we will concentrate on the late quintessence 
model. 

Even though we use the full set of the above equations for our calibration, 
we investigate the eflFects of coupling on the matter power spectrum with 
some approximation in what follows. It is well known that averaging out the 
small and oscillatory 5Q" and 6Q' is a good approximation for all scales [9] . 
Prom this fact, we obtain the approximate expression of 5Q from Eq. (6) 

5Q ~ r , — , 8 

where niQ = d^V/dQ'^. We obtain the evolution equation of 5c from Eqs. (4) 
and (5) by using Eq. (8), 

5'^ + h[i + nc^J^{l + 0JQ)^Q] 6', - ^n^nc{l + 2nl) 5e ^ , (9) 
where Q,q is the quintessence energy density relative to the critical density 



and we use Q' = 'H^?,{1 + ujq)Q.q [4]. The evolution of the linear pertur- 
bation of the baryon ^5 which is not coupled to the scalar field is the same 
as the above equation (9) except that now the coupling terms are absent: 

5'{, + nS'^ - ^-H^^cSc ^ . (10) 
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It is convenient to rewrite the above equations (9) and (10) in terms of 
a; = In a as 



+ 



d6c 
dx 



dx'^ 



+ 



d6b 
dx 



31 
2V 



where Qq ~ 1 — ilc — 



1+ ' 



6-0, 



(11) 
(12) 

(13) 



In Eq. (12), we define a baryon bias factor b hy 5b = b6c with a smah time 
variation {i.e. we ignore db/dx and d'^b/dx'^ terms). If we further use the 
assumption that the hnear growth of the CDM is given by 6c oc e™^ with 
m < 1, then we will obtain the analytic form of b 



l + 2ni 



2ncmy3(l + ujQ)il.Q 
30^ 



(14) 



Here we use the ansatz for the coupling between CDM and DE as e"""^ . In 
our model, the scalar field evolves from —5 to during 10"^'^ < o < 1- We 
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Figure 2: a) Evolutions of 6^ in the ACDM, QCDM, and cQCDM models 
(from top to bottom), b) Evolutions of the equation of state of Q field, loq, 
for the coupled and noncoupled cases (from top to bottom). 



put a limit on the magnitude of the coupling constant \nc\ < 0.01, where 
the choice of the sign of ric is still arbitrary. The effective mass of the 
CDM varies at most around 5 %. rric {equally, pc) increases (decreases) as it 
evolves to the present for the positive (negative) ric- Thus, the evolutions of 
the background quantities are slightly changed dependent on the sign of Uc- 
The effect of the coupling on the h also depends on the sign of ric as given 
in Eq. (14), which is used to estimate h as shown in Fig. 1. In the left panel 
of Fig. 1, we show the evolution of h when ric = +0.01. For the positive 
the last term in Eq. (14) is negative and its magnitude is bigger than 2n^ for 
the given model. Thus, 5 > 1. The evolution of h for ric = —0.01 is depicted 
in the right panel of Fig. 1. In this case, the last term in Eq. (14) is positive 
and b is always smaller than 1. Thus, we are able to constrain not only the 
magnitude but also the form of the coupling between CDM and DE from 
accurate observations of the baryon power spectrum. A similar but slightly 
different conclusion was drawn in Ref. [10]; however, their conclusion is only 
true for the tracking region solutions. 



First, we study the CDM density fluctuation 5c given in Eq. (11) for dif- 
ferent cases. We denote respectively the cosmological model including CDM 
component with the cosmological constant as ACDM, with the non-coupled 
Q field {ric = 0) as QCDM, and with the coupled Q field as cQCDM. We also 
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denote 5c for each model as 5^^^^, 5^^^^, and df^^^. 5^^^^ has the 
same evolution equations as those of the ACDM model except the difference 
in Q,Q given in Eq. (13). ujq in the QCDM model changes from 1/3 in the 
radiation-dominated epoch (the so-called "early tracking region") to around 
— 1 at present. Thus, > oj\ = ~1 during the entire epoch. This causes 
the suppression of 6c in QCDM model compared to ACDM model when we 
use the same present values of the cosmological parameters. We illustrate 
this in the left panel of Fig. 2. The diamond and the rectangular points cor- 
respond to 6c'" and 6^'-''^^^ , respectively. We also compare 6'^'~^^^ and 
^cQCDM _ -^Vhen the CDM is coupled to Q field, the scaling of is changed 
as given in Eq. (2). Also lcq is increased during the matter-domination 
epoch in the cQCDM model as shown in the right panel of Fig. 2. The dot- 
dashed and the dashed lines depict ujq when Uc = 0.01 and 0, respectively. 
Thus, this causes slightly further suppression of Sc in the cQCDM model. 
However, if we want to compare (JcQCDM ^h]^ ^QCDM ^j^g relevant sub- 
horizon scale, then we should constrain their evolutions at late times a > 0.1 
{equally, x > —2.3). But then ujqs of the two models are almost identical 
as shown in the right panel of Fig. 2 and the discrepancy between (J^QCDM 
and (jQC'^^ becomes negligible. If we use the definition of the baryon bias 
factor 6=1^, then we would obtain that bcQCDM > bqcDM — b\CDM — 1 
for the same present values of the cosmological parameters. 

Now we consider the previously mentioned differences in the two-point 
correlation functions in different models. Instead of using the above approx- 
imations, we will run the numerical evolution of the full system of equa- 
tions. The coupling of the quintessence to the dark matter modifies both 
the turnover scale and the amplitude of the matter power spectrum. How- 
ever, the shift in the turnover scale and the suppression in the amplitude 
of the matter power spectrum (defined by Pik) = = (27r^)//c^A^(A;)) 

due to this coupling are too small when we limit the coupling strength 
to be consistent with observations. Thus, this gives us the motivation to 
probe the coupling effects on the baryon acoustic peak in the correlation 
function, ^{s) = J A'^(k)jQ{ks)dlnk, where Jq is the spherical Bessel func- 
tion. In order to avoid the damping effect of the BAO signature in the 
matter power spectrum due to nonlinear clustering, we put the limitation 
of A; < 0.2/iMpc~^ in our correlation function calculation. We show the 
correlation function times the comoving separation square (s^) in Fig. 3. 
The solid, dashed, and dotted lines correspond to the ACDM, QCDM, 
and cQCDM models, respectively. The cosmological parameters that we 
use in this figure are Hq = 71 km/sec/Mpc, Q.^, = 0.047, Qc = 0.211, 
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Figure 3: Matter correlation functions for ACDM, QCDM, and cQCDM 
models. 

= 0.742, and the galaxy bias factor bgai = 1.9. We also normalize 
the matter power spectrum to cg = 0.788 to be consistent with the WMAP 
and the Luminous Red Galaxies (LRG) observations [11]. Note that the in- 
fluence of DE on the present value of dg = \/ Jq W^{kR8)A'^(k)dlnk, where 
W{x) = 3(sinx/x^ — cosx/x^) is the Fourier transform of a top-hat window 
function and Rs = 8/i~^Mpc, is well studied [12, 13, 14, 15]. 

In Fig. 3, we are able to clearly see the differences of the correlation 
functions for the different models. The first peak in the correlation func- 
tion corresponds to the turnover scale which is related to the scale factor 

~ 1/3500 when the radiation and the matter densities are equal. The 
ACDM and QCDM models have the same agq for the same set of cosmo- 
logical parameters while their comoving distances of the matter and radi- 
ation equality, Sgg, are 19.17 Mpc/h and 19.60 Mpc//i, respectively. The 
discrepancy comes from the fact that the comoving distance to the a^g is 

/o"' ^^Fy ~ Iz^q H(z') ^^'^ models have slightly different H{z). Also, 

the location of Seq in the cQCDM model is different from that in the QCDM 
model due to the change in the scaling of the CDM density. In our model, 
e > and it causes the delay of the radiation and matter equality epoch 



QCDM 
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(Og^'"^^ > a'^q^'^^). Thus, the comoving distance s^q is shifted to about 
20.87 Mpc//i in the coupled case. 

We have already shown for the relevant A;-modes, > ^QCDM ^ 

gcQCDM £qj. same present values of the cosmological parameters. Thus, if 
we normalize the power spectra of the latter two models on cluster scales erg, 
they acquire a larger amplitude of primordial fluctuations compared to the 
ACDM model. The BAO bump in the correlation function of both QCDM 
and cQCDM models are larger than that of the ACDM model. However, 
the reason for the enhancements of the BAO bumps in both models are 
different. The enhancement in the QCDM model compared to the ACDM 
one is due to the choice of the larger amplitude of primordial fluctuations. 
The enhancement in the cQCDM model is due to the coupling between Q 
and CDM. We observe this effect in Fig. 3. If we compare the BAO peak of 
cQCDM with that of ACDM or QCDM, then we observe that it is enhanced 
in the cQCDM model. This is consistent with our early explanation that the 
baryon bias factor in the cQCDM model is enhanced as given in Eq. (14). 
The overall amplitudes of all three models in Fig. 3 are smaller than those 
given in Ref. [16]. This is due to the different choices of the galaxy bias 
factor. We use the LRG galaxy bias factor bgai = 1.9 [11] whereas the au- 
thors in Ref. [16] claim that they use the scale dependence bias factor which 
is however not given therein. However, both the shapes of the correlation 
functions and the amplitudes of the BAO bumps of the QCDM and cQCDM 
models fit to the data points better than the ACDM model. As we men- 
tioned above, if we choose a negative n^, then the amplitude of the BAO 
bump in the cQCDM model is decreased due to anti-biasing of baryons 
with 6 < 1. This effect mimics the nonlinear effect [17]. It will suppress 
the amplitude and shift the location of the BAO bump, and thus making 
harder for us to fit the data. However, we limit the calculation in the linear 
regime and this effect is irrelevant for our consideration. In the Figure 3 
of Ref. [16], it is claimed that ^ch^ = 0.13 shows the better fit to the data 
compared to the case using Qch^ = 0.12. However, the amplitude of the 
BAO bump with Qch'^ = 0.13 shows the bigger discrepancy with the data 
than that with ^ch^ = 0.12. We are able to give a better direction for fitting 
both the shape and the amplitude of the correlation function in the QCDM 
and cQCDM models. We have done a simple analysis and found that 
the values arc 21.6, 21.7, and 23.4 for the ACDM, QCDM, and cQCDM 
models, respectively. Although the improvement has low statistical signifi- 
cance, the species-dependent coupling effects to the correlation function in 
the coupled quintessence model may be potentially important and are being 
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further studied. Note that there is a slight shift in the location of the BAO 
peak in the QCDM model. This is due to the change in H{z) in the QCDM 
model compared to the ACDM model. This affects the sound horizon that 
is given by Is = -rff-rdz, where the sound speed of the baryon-photon 

fluid Cg = l/-\/3(l + 3pij/4pr) is the same for every model. Hence, for each 
model we have a different H(z) and Ig will vary. However, the difference is 
quite small. 
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